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1 = 2 Pion Scattering Length from Two-Pion Wave Function * 
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a Center for Computational Physics, University of Tsukuba,Tsukuba, Ibaraki 305-8577, Japan 
b Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 

We present a preliminary report on a calculation of scattering length for 1 = 2 S-wave two-pion system directly 
from the two-pion wave function. Results are compared with those calculated from the time dependence of two- 
pion four-point functions. Calculations are made with an RG-improved action for gluons and improved Wilson 
action for quarks at a' 1 = 1.207(12) GeV on 20 3 x 48 and 24 3 x 48 lattices. 



1. Introduction 

By now the standard procedure employed for 
calculating the scattering length of hadrons from 
Euclidean lattice QCD is to use the finite-size 
method proposed by Liischer, which relates the 
finite volume shift of energy eigenvalues AE to 
the scattering lengths. The formula reads [7] 



ao/(Lir) — —x — A ■ 



B 



(i) 



where x = AE ■ 2m 7r L 2 /(47r) 2 , and A = 
—8.9136 . . . and B = 95.985 . . . are geometrical 
constants. The energy shift AE is extracted from 
the time dependence of two-pion four-point func- 
tions. In this way a number of calculations has 
been made for the 1 — 2 S-w&ve two-pion scat- 
tering length using the staggered [1,2], the stan- 
dard [2-5] and the improved Wilson fermion ac- 
tions [6]. 

It should be noted that the condition R < L/2 
is assumed for the two-pion interaction range R 
and the lattice volume L 3 in derivation of the 
formula (1) [7]. So far there have been studies of 
the lattice volume dependence of the scattering 
length, but no direct investigation of the interac- 
tion range R. It is extremely important that we 
examine the validity of the necessary condition 
for the Liischer's formula in our current lattice 
simulations. This can be done by investigating 
the two-pion wave functions, which is one of the 
purpose of this article. Similar work has been 
carried for 2 dimensional XY model [81. 



Once one has access to the wave functions, 
one can try to extract the scattering length from 
them. This is the second, and perhaps more in- 
teresting, purpose of this article. 

2. Method 

Our idea is based on the derivation of the 
Liischer's formula (1). Liischer found that the 
two-pion wave function <p(x) on a finite pe- 
riodic box L 3 satisfies the following effective 
Schrodinger equation [7]. 



(A + k 2 ) ( f ) (x) = / d 3 y U k (x,y)<j>(y) 



(2) 



where x is the relative coordinate of the two pion 
and k 2 is related to the two-pion energy eigen- 
value through E = 2 • y / m 2 r + k 2 . The function 
Uk (x, y) is the Fourier transform of the modified 
Bcthe-Salpeter kernel introduced in [7]. It is non- 
local and generally depends on energy. 

In the derivation of the formula it is assumed 
that Uk(x,y) ~ for R < \x\, \y\. In this region 
the wave function satisfies the Helmholtz equa- 
tion (A + k 2 )4>(x) — 0. The general solution of 
Helmholtz equation on a finite periodic box L 3 
can be written as 



</>(£) = £ e 1 - 

n£Z 3 



i/N 



/{n 2 -k 2 -N 2 ) 
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up to over all constant, where N = L/(2n). As- 
suming the interaction range is smaller than one 
half of box size, ie. R < L/2, and expand- 
ing the general solution (3) in terms of spheri- 
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Figure 1. Two-pion wave function on 20 3 lattice 
at t = 30 and ml = 0.27GcV 2 . 
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Figure 2. A<j>(x,t)/<j>(x,t) ( k 2 ) on 20 3 lattice 

at ml = 0.27GeV 2 for several t. 



cal Bessel ji{x) and Noeman ni(x) functions for 
R < \x\ < L/2, we obtain 

(j>{x) = a (k) ■ j (k\x\) + Po(k) ■ n (k\x\) + • • • (4) 

where neglected terms are contributions from 
states with angular momentum I > 4. The expan- 
sion coefficients ao(k) and (3o(k) yield the scatter- 
ing phase shift in infinite volume by tan^o(fc) = 
Po(k)/ao(k). In particular for the lowest en- 
ergy state of the two-pion, it gives the scattering 
length by o (k)/a o (k) = a -k + 0(k 2 ). The con- 
stant Po(k) are also related to ao(k) geometrically 
and the relation leads to the Liischer's formula 
(1). 

In this work we define the two-pion wave func- 
tion by 

<f>(S,t) = ^2(ir{R{g)+Jl,t)ir{Jl,t) S) (5) 

where R is an element of cubic group, and sum- 
mation over R and X projects out the A + sector 
of the cubic group and that of the zero center of 
mass momentum. In order to enhance signals we 
use a source constructed with two wall sources 
given by S = W(*i)W(t 2 ), with W(t) the wall 
source at t, in Coulomb gauge fixed configura- 
tions. 



3. Results 

We work in quenched lattice QCD employing 
an RG-improved action for gluons at /? = 2.334 
and an improved Wilson action for quarks at 
Csw = 1-398. The corresponding lattice cutoff is 
estimated as 1/a = 1.207(12) GeV from m p . The 
volumes of lattices (number of configuration) are 
20 3 x48 (150) and 24 3 x48 (100) which correspond 
to 3.26fm 3 and 3.92fm 3 in physical units. Quark 
masses are chosen to be m^GeV 2 ) = 0.27, 0.44, 
and 0.74. Quark propagators are solved with the 
Dirichlct boundary condition imposed in the time 
direction and the periodic boundary condition in 
the space directions. The sources are set at t\ = 7 
and ti = 8 for both lattices. 

In Figure 1 we show the ratio of the wave func- 
tion (p(x, t)/4>(xo, t) on a 20 3 lattice at t = 30 and 
m 2 = 0.27GcV 2 . Here we choose |x = 7.6 for 
normalization point. We find that the statisti- 
cal error is very small. Figure 2 shows the ratio 
A<j)(x,t)/(j)(x,t) at the same parameters for sev- 
eral t. This ratio should take a negative constant 
— k 2 for the 1 = 2 two-pion system for R < \x\. 
Clear plateau appears for \x\ > 7.5 at t > 30 as 
shown in the figure. This indicates that the neces- 
sary size for application of the Liischer's formula 
is L > 2 x 7.5a (= 2.45fm ). 

We try to extract the value of k 2 from the 
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Figure 3. 1 = 2 two-pion scattering length 
an/m^ (1/GeV 2 ) obtained from the time corre- 
lation functions and the wave functions. 



wave function by fitting the data with the func- 
tion defined by (3). The results are plotted by 
large cross symbols in Figure 1, where we choose 
\x\ = 7.8 — 16 as the fitting range. The fit works 
very well. The result for k 2 is 2.67(12) x 1(T 3 
GeV 2 , which is consistent with that obtained 
from the time dependence of the two-pion four- 
point function, 2.62(23) x 1(T 3 GeV 2 . We check 
the consistency between the two determinations 
at all simulation points on 20 3 lattice in Fig- 
ure 3, where the scattering length is calculated 
by substituting k 2 to the Luscher's formula (1). 
We observe that the statistical error of the scat- 
tering length obtained from the wave function is 
smaller than that from the time correlation func- 
tion. This feature was also found for the 2 dimen- 
sional XY model [8]. 

The wave function on a larger lattice of 24 3 at 
to 2 = 0.74GeV 2 and the results of the fitting are 
shown in Figure 4, where we choose t = 30. We 
find a systematic discrepancy between the data 
and the fit for large |x| region (|a?| > 15). Fur- 
ther there is a large discrepancy in the scattering 
length obtained from two methods as shown in 
Figure 3. 

A possible cause of the discrepancy is smallness 
of the time extent. The strength of the two-pion 
interaction on a finite L 3 periodic box is propor- 
tional to 1/L 3 . If we consider a time-dependent 
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Figure 4. Two-pion wave function on 24 3 lattice 
at t = 30 and ml = 0.74GeV 2 . 



effective Shrodinger equation, the wave function 
satisfies a heat diffusion equation. We expect that 
a larger time extent is necessary to obtain the cor- 
rect wave function for a larger volume. In order 
to confirm this anticipation, we are going to ex- 
plore simulations on a temporal extent T = 80 in 
our next study. 
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